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ines won Ob SOll CONSOLIDATION 
AND TESTING OF FOUNDATION SOILS 


BY THE DIVISION OF TESTS, BUREAU OF PUBLIC ROADS 


Reported by L. A. PALMER, Associate Chemist, and E. S. BARBER, Junior Highway Engineer 


pression test described in this report were con- 

structed by the Bureau of Public Roads in 1926 in 
accordance with plans furnished by Dr. Charles 
Terzaghi. Investigation of testing technique and 
practical utilization of the resulting data have been 
included in the Bureau’s program of research on soils 
for highway purposes continuously since that time. 

As part of this program, the Bureau cooperated with 
the Massachusetts Institute of Technology in research 
on soil mechanics from 1926 until the fall of 1932. Dr. 
Terzaghi was retained by the Bureau as research con- 
sultant on the work from the fall of 1926 to the fall of 
1929. Dr. Arthur Casagrande was employed by the 
Bureau as research assistant and was stationed at the 
Massachusetts Institute of Technology throughout the 
duration of this cooperative research work. 

Other parts of the Bureau’s general program were 


Sica: DEVICES used in performing the com- 


carried on at the Bureau’s laboratory at Arlington, Va., | 


and through cooperation with various State highway 
departments and Federal bureaus. 


Special studies have been made at Arlington to 
determine the magnitude of experimental errors result- 
ing from variables in testing apparatus and procedure, 
and to determine the agreement between deformations 
observed in tests with those computed in accordance 
ene theory of consolidation on which the tests are 
based. 


The compression test and theory of soil consolidation 
were described by Terzaghi in 1927.1 Two other 
reports of distinctive value to those concerned with the 
practical utilization of the compression test data have 
also been published.’ 


On the basis of the Bureau’s investigations of the test 
itself, experience in practical use of the test data, and 
reports just referred to, it is now possible (1) to state 
definitely that the apparatus and testing procedure 
described briefly in this report have been found satis- 
factory, and (2) to present a practical method of esti- 
mating that part of the total settlement of soil caused 
by the loss of water forced vertically out of saturated 
compressible soil strata in certain soil profiles. 

Settlements of soils in profiles differing from those 
specifically described, or produced by plastic flow, flow 
of water laterally, or causes not discussed, are beyond 
the scope of this presentation. 


COMPRESSION DEVICE VALUABLE AID IN STUDYING SOIL 
COMPRESSION 


The theory of consolidation of earth materials in- 
volves two problems: (1) Determination of the distri- 
bution of stresses; and (2) computation of the displace- 
ments caused by these stresses. In dealing with many 
engineering materials and structures the essential 


1 Principles of Final Soil Classification, PUBLIC ROADs, Vol. 8, no. 3, May 1927. 

2 Report of the Special Committee on Earths and Foundations, Proceedings of 
the American Society of Civil Engineers, vol. 59, no. 5, May 1933, and discussion by 
William P. Kimball, Proceedings of the American Society of Civil Engineers, vol. 
59, no. 6, August 1933. 
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problem is to estimate the stresses, since nearly all 
structural materials possess, within limits, the property 
of elasticity. With soils, however, the problem is 
quite complicated for they have both plastic and elastic 
properties; the materials are essentially heterogeneous; 
and water is usually present. If the water is free to 
flow from the compressed soil most of the settlement 
results directly from the reduction in voids between 
solid particles, an occurrence that, in this paper, is 
considered apart from elastic and plastic deformations. 

One of the most useful tools in the study of soil con- 
solidation is the laboratory compression device designed 
by Dr. Terzaghi. The compression tests indicate 
within a comparatively short time in the laboratory 
the consolidation characteristics of soils upon which 
depend the rates of settlement and total settlements 
of foundations. Thus, compression tests supplant 
guesswork and “rule of thumb” methods with reliable 
information for use in the design of foundations. 

Figure 1 illustrates the essential features of the test- 
ing device used by the Bureau of Public Roads. An 
undisturbed sample of soil with natural structure and 
moisture content, as nearly as can be obtained, is 
carefully placed between two porous stones, which 
act as filters during the test. Pressures are then ap- 
plied, ranging from zero to the maximums that are 
expected in service. A burette connected with the 
apparatus as shown in figure 1 affords a means of de- 
termining the permeability of the sample at any time 
during the test. 

The sample during test is representative of a stratum 
of saturated, compressible soil in the earth’s crust, 
sandwiched between two layers of sand or other rela- 
tively permeable material. 

The total amount of consolidation and the speed at 
which it takes place are controlled by: (1) The moisture 
content of the soil; (2) the velocity with which water 
flows vertically through the pores of the sample to the 
filters above and below in the laboratory, or through 
the compressible layer to the sand strata in the field; 
and (3) the frictional resistance of the soil particles to 
consolidation. 

At the instant of load application, the water in the 
compressible sample or layer is considered as furnishing 
all of the resistance to consolidation. Thus, if a suitable 
pressure gauge were connected with the sample in the 
laboratory, it would show that, at the instant of load- 
ing, the water pressure in the sample is practically 
equal to the pressure applied to the sample. It was 
stated during the first conference on soil mechanics and 
foundation engineering, held at Harvard University in 
June, 1936, that experiments in Sweden disclosed that 
the water pressures in compressible understrata were 
approximately equal to the pressures produced by the 
weight of superimposed fill materials just at the time 
construction of the fill was completed. 

As the pressure on the soil continues, water escapes 
from the sample or field stratum, causing the pressure 
to be gradually transferred to the soil particles or 
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Figure 1.—Dervick For DETERMINING THE CONSOLIDATION 
CHARACTERISTICS OF SOILS IN THE LABORATORY. 
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FiguRE 2.—PERMEABILITY, COMPRESSION, AND CONSOLIDATION 
CURVES FOR A SAMPLE OF COMPRESSIBLE SOIL. 


“soil skeleton.” Eventually, equilibrium is reached 
and further consolidation ceases when the water pres- 
sure becomes equal to zero and the soil skeleton pro- 
vides all of the resistance to consolidation. 

The sample in the compression device is subjected to 
pressure that is increased by increments. Each in- 
crement of pressure acts until apparent equilibrium is 
attained before the next increment is added. Water 
escapes under pressure, entering the porous stones and 
flowing from the outlets as indicated in figure 1. 

Consolidation characteristics and permeability of the 
laboratory specimen are recorded and the data are 
plotted as shown in figure 2. The load-compression 
curve shows the voids ratios at which the soil sample 
attained apparent equilibrium under different pressures. 
The voids ratio is defined as the ratio (by volume) 
of voids to solids for a given sample of soil. 

Apparent equilibrium in the sample is reached in a 
relatively short time, and for all practical purposes is 
true equilibrium. Actually, however, the attainment 
of true equilibrium would require infinite time. For 
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values of the pressure, p, greater than about one-tenth 
kilogram per square centimeter, the load-compression 
curve, as shown, is expressed by the equation,® 


ée=1.050—0.442 log ps. - 1 (1) 


where e is the voids ratio. The permeability-voids 
ratio curve (fig. 2, upper) shows the coefficient of 
permeability of the sample at different stages of 
compression. 


COEFFICIENTS OF PERMEABILITY AND COMPRESSIBILITY 
DETERMINED 


The coefficient of permeability is defined as the 
velocity of flow when the hydraulic gradient is unity, 
that is, when the drop in hydrostatic pressure in grams 
per square centimeter between two points in a soil 
column equals the distance in centimeters between the 
points. The usual practice is to record and plot the 
coefficient of permeability in centimeters per second, 
asin figure 2. However, for the purpose of simplifying 
subsequent computations the coefficient of permea- 
bility, as used in this report, is expressed in centimeters 
per minute. The average time-consolidation curve 
(fig. 2, lower) shows the rate of consolidation under 
any given pressure. ‘This curve is the average of the 
individual time-consolidation curves for different load 
increments. Essential features of the test procedure, 
the equation for determining the coefficient of perme- 
ability, and the use of average time-consolidation 
curves have been discussed previously.* 

Both the consolidation characteristics of the soil 
strata beneath a given footing and the vertical pressures 
at various depths produced by the superimposed load 
must be known. In order to solve the many different 
problems involved in computing earth pressures, one 
must be familiar with the mathematical theory appli- 
cable to all types of problems requiring the calculation 
of earth stresses and their distribution. The treatment 
of special cases is then simplified. 

The coefficient of permeability, k, is the velocity of 
flow at a hydraulic gradient of unity. If the thickness 
of the soil sample, with a voids ratio of e, is d, and if the 
voids ratio could be reduced to zero, the thickness 
would become dp, that is, 


d=d, (7 ¢) .-82.22 eee (2) 


The rate at which water would permeate the soil 
sample if the difference between the two water levels 
were do, rather than d, is designated as ky. It has been 
shown ' that the rate of permeation increases as the 
hydraulic gradient increases so that from this fact and 
equation (2), 


kod 
kein lh¢ 2 2s eee (3) 
that is, 
k 
es pay, ala (4) 


1 Principles of Final Soil Classification, by Dr. Charles Terzaghi, PuBLic Roaps, 
vol. 8, no. 3, May 1927. 

§ When plotted to a semilogarithmic scale, e being plotted to the natural scale, the 
load-compression curve is generally a straight line for loads greater than 0.1 kilogram 
per square centimeter. The general equation for the straight-line portion of the 
curve is 


e=B—Z log p 


where B is the intercept of the straight-line portion of the curve at the 1.0-kilogram- 
per-square-centimeter ordinate and — Z is the slope of the semilog plot of the curve, or 
range in value of e for unit change of log p. 

‘A Method of Predicting Settlement of Fills Placed on Muck Beds, by F. A. 
Robeson, PUBLIC Roaps, vol. 16, no. 12, February 1936. 
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It is convenient to take e as the average voids ratio of 
a soil sample compressed from an initial voids ratio, 
é;, to an equilibrium value e:, so that 


On substitution of this value for e in equation (4), 
one obtains 
2k 


3 fe her 6) 


The coefficient of compressibility, a, is the average 
decrease in voids ratio per unit increase of pressure. 
If the initial voids ratio, ¢;, is the condition at equilibrium 
under a pressure of p, grams per square centimeter and 
€: is the corresponding condition at a pressure p2 grams 
per square centimeter, then 


ko 


The values of both a and k decrease with the voids 
ratio, but at different rates. However, the change in 


the value of the ratio,“ is fairly small. The ratio, = 


is called the coefficient of consolidation, ¢. From 
equations (6) and (7) there results from substitution for 
ky and a, 





hee PsP) 2k 
peeks vrei Ciena \ f. (8) 


a €;— 62 

The utility of this equation may be illustrated by use 
of the data plotted in the upper curves of figure 2. 

From this curve, the pressures, p,; and p2, correspond- 
ing to ¢;=1.00 and e.=—0.80, are 1,300 and 3,680 grams 
per square centimeter, respectively. The average voids 
ratio is 0.90 and from the permeability curve, the 
coefficient of permeability corresponding to a voids 
ratio of 0.90 is 4.6 107-8 centimeters per second. 


__3,680—1,300., 9.2107 


Then c=~7 00—0.80 X3+1.00-0.80 


centimeter per minute. 

Similar computations can be made, taking different 
values for e, and e, and the corresponding pressures, 
p, and p, from the curve, finding the coefficient of 
permeability corresponding to the average voids ratio 
as illustrated above, and substituting these values in 
equation (8).° By this procedure, the values of c given 
in table 1 are computed. These values are plotted 
against the corresponding voids ratios in figure 3. 


X60=0.0173 


5 The pressures for voids ratios of 0.6, 1.2, and 1.4 were computed from the relation, 
€=1.050—0.442 log p:p. 
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Figure 3.—RELATIONS BETWEEN COEFFICIENT OF CONSOLIDA- 
TION AND Vorps Ratio. 


TABLE 1.—Computed values of the coefficient of consolidation, c 


























- - Pressures corresponding | Coefficient of 
Woids favics to voids ratios permeability, k,! Qoemcient 
corresponding of consolida- 
2 ” tion, ¢c 
171762 
2 
é1 €2 pi p2 
Grams per | Grams per ; 
square centi- | square centi- | Centimeters per | Centimeters 
meter meter second per minute 
1, 40 1, 20 161 458 22. 61078 0. 0088 
1. 20 1.00 458 1, 300 10, 21078 . 0123 
1.00 . 80 1, 300 3, 680 4.61078 . 0173 
. 80 - 60 3, 680 10, 400 1. 610-8 . 0190 


This paper introduces only the most essential steps 
in some of the mathematical derivations used both for 
the computation of earth pressures and for computation 
of formulas applicable to the consolidation of relatively 
impervious, saturated strata caused by these pressures. 
For readers more concerned with direct application of 
the fundamental principles and formulas than with 
their development, a discussion of the use of the theory 
is presented in detail in part 1 of the paper. The mathe- 
matical derivations are included as part 2 for those 
interested also in the theoretical development. 
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PART 1.—APPLICATION OF THEORY OF SOIL CONSOLIDATION 


The manner of applying the test data in practice 
varies, depending upon the pressure distribution condi- 
tions beneath the loaded area, the thickness of the 
compressible layer, and whether it is free to drain both 
at top and bottom, as in the compression test, or from 
only one of the bounding surfaces. 

When a Joad is applied on the top of a compressible 
soil stratum, conditions may be such that the pressure 
resulting from the load is the same from top to bottom 
of the stratum, that is, the pressure distribution is 
uniform. On the other hand, the conditions may be 
such that the pressure varies, 1. e., increases or decreases 
as the depth increases. In this paper the pressure 
change per unit of depth is considered to be either zero 
or constant. If it is zero, the pressure distribution js 
uniform; if it is constant, the pressure distribution is 
either triangular or trapezoidal. The cases of pressure 
distribution are illustrated in figure 4 and are as follows: 

Case 1.—Rectangular pressure distribution. The 
pressure is uniform from top to bottom of the compres- 
sible stratum. 

Case 2.—Triangular pressure distribution. The pres- 
sure varies from zero at the upper bounding surface of 
the compressible stratum to a finite value greater than 
zero at the lower bounding surface. 

Case 3.—Triangular pressure distribution. The pres- 
sure varies from a finite value greater than zero at the 
top of the compressible layer t6 a pressure of zero at the 
bottom. 

Case 4.—Trapezoidal pressure distribution. The 
pressure increases with increasing depth in the com- 
pressible stratum. The pressures at both top and 
bottom surfaces are finite values, each greater than zero. 

Case 5.—Trapezoidal pressure distribution. The 
pressure diminishes with increasing depth in the com- 
pressible stratum. The pressures at both top and 
bottom surfaces are finite values, each greater than 
zero. 


FORMULAS DEVELOPED FOR VARIOUS TYPES OF LOADED AREAS 


The pressure distribution throughout the compres- 
sible layer depends on the shape and magnitude of the 
loaded area at the ground surface and on the load 
distribution over this area. In the examples cited in 
this paper the load distribution is taken as uniform 
over the loaded area. 

The five types of loaded areas considered are: 


A.—A loaded area of great length and width. 

B.—A spot load. 

C.—A relatively short and narrow, rectangular 
loaded area. 

D.—A relatively narrow loaded area of great 
length. 

E.—A circular loaded area. 


The pressure per unit area at any given depth below 
a loaded area of great extent is taken as equal to the 
unit pressure on the loaded surface plus the weight of 
superimposed earth above the point in question. No 
special formula is required for computations in this 
simple case. 

To compute the pressure below a spot load, use is 
made of the equation, 











JM 





CASE ] CASERzZ CASE 3 CASE 4 CASE 5 
Figure 4.—Five Casts or PRESSURE DISTRIBUTION IN Soin 
STRATA. 

where 
(10) 


P=the spot load; and 


pz=the vertical pressure per unit 
area at a vertical distance, 2, 
below a spot load on the sur- 
face and at a horizontal dis- 
tance, 7, from the load. 

For any distance on the vertical axis below the load, 

: Le 
ris zero and Le oa 

Equations (9) and (10) can be used to determine the 
vertical pressure produced by any load, uniformly 
distributed over a rectangular area, provided the ratio 
of the longer side of the rectangle to z, the vertical 
depth to the point at which p, is to be computed, is 
not greater than one-half. This condition is illustrated 
in ficure 5. If, however, this ratio exceeds one-half, 
then the procedure is to subdivide the rectangle into 
rectangular elements (all of equal area) such that the 
ratio of the longer side of a single element to the ver- 
tical distance to the point in question is not greater 
than one-half. If the point at which it is desired to 
know the vertical pressure is at a very considerable 
depth below the loaded footing, then the procedure of 
subdivision may be unnecessary even in the case of 
large loaded areas. For example, the base of the 
Washington Monument is 125 feet square. At any 
point 250 feet beneath the base the ratio is 5 
and hence the entire loaded area may be considered as 
a point load at this depth. 

When subdivision is necessary it is considered that 
the load on each element is concentrated at its center. 
To determine the unit pressure on the vertical center 
axis, the horizontal distance r from the center of each 
element to the center of the footing is computed and 


the ratio, % is evaluated. This is substituted in the 


expression for K, equation (10), and K is computed for 
the element of area. Having done this with all of the 
different elements of area, the AK values are added 
together and the sum is substituted in equation (9), the 
value for P in this equation being the total load on each 
individual subdivision of area. This procedure will be 
further illustrated later. 

For a relatively narrow loaded area of great length, 
take a vertical cross section of the loaded strip. At a 
point on the vertical center axis, p, is computed from 
the equation, 
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Figure 5.—REcTaNGuLAR AREA LOADED UNIFORMLY. (AREA 
or Typp C, ILLUSTRATING THE PRINCIPLES oF BotH TypEs 
B anv C.) 


p:="(a+sin 5h 1 a Mac a 9 pa eee pea 


a= twice the angle, in radians, whose tangent 
is half the width of the strip divided by 
the vertical distance z to the point in 
question (see fig. 6); a p=the load 
per unit area of the footing. 


where 


To compute the maximum vertical pressure, p., at 
any point below a circular footing and on the center 
line (see fig. 6), use the equation, 


eee eG py eee ee Pe 2. fe (12) 


where p=the load per unit area of the footing 


and 8=the angle whose tangent is (see fig. 6). 


RATE OF SETTLEMENT COMPUTED BY SIMPLE FORMULAS 


If the compressible layer is sandwiched between two 
permeable sand-gravel layers, one above and one below, 
water can escape in two directions. If the compres- 
sible layer is overlaid with sand and rests on rock, 

water can escape from but one face of the stratum. 
There is a double filter in the first case and a single 
filter in the other. In either case, under a given 
pressure the thickness of the compressible layer is 
reduced to the same extent, ultimately. However, the 
rates of reduction in thickness of the compressible 
layer will obviously be different in the two cases. Con- 
sequently, the settlement during a relatively short time, 
say in 6 months or a year, will be greater if filtration 
can occur at both the upper and lower boundaries of 
the compressible stratum. The filtration effects a 
diminution in thickness and the earth’s surface above 
the compressible stratum settles. 

The total amount of vertical consolidation Q is 
computed by substitution of the data shown by the 
load-compression curve in the following expressions: 


For soil sample, 


Y= besa: a a gee ak Te 


for soil stratum, 
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where e,=voids ie of the sample, taken as the aver- 

age voids ratio in the compressible 
stratum prior to the application of addi- 
tional load at the surface; 

€2—voids ratio of the sample taken as the average 
voids ratio in the compressible str atum 
upon reaching equilibrium after the appli- 
cation of additional load at the surface; 

d,=thickness of soil sample at apparent equilib- 
rium at voids ratio, ¢,; and 

D,=the thickness of the soil stratum at voids 
ratio, é1. 


The rates at which soil strata settle are computed 
from test data and field observations by simple formulas. 

The theoretical time-settlement relation for soil 
samples or strata sandwiched between two porous 
layers and with either rectangular, triangular, or 
trapezoidal pressure distribution or for strata with 
rectangular pressure distribution and with a porous 
layer only at the upper face is oe by the expression, 


n=1—S(¢ ee se We ag dere (15) 


where q;=percentage of ultimate consolidation at the 
end of a definite period of time, ft, of load 

application; 

e= Naperian base. 

The exponent, N, is related to the thickness 
of the soil sample, the period of load appli- 
cation, f, and the speed of egress of water 
from the soil sample or stratum. 
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Experience has shown that the time-consolidation 
curves furnished by the compression test on compressi- 
ble soils (see fig. 2, lower curve) conform in shape with 
the theoretical curve of equation (15) from g, equals 
zero to gq, equals 90 percent on the theoretical curve. 
The experimental values of q attributable to per- 
meability alone, are usually from 60 to 100 percent of 
the actual ultimate consolidation as obtained in the 
laboratory with the given loading conditions. Only 
laboratory values of g, within the range in which the 
experimental curve agrees with the curve showing the 
theoretical relationship are valid for computations 
involving the use of equation (15). Beyond the limit 
of consolidation as effected by the egress of water, 
further consolidation of most soils is caused by particle 
rearrangement and does not result from permeability 
as assumed in the derivation of equation (15). 


TIME-CONSOLIDATION RELATIONS GIVEN FOR UNIFORM, TRIANGU- 
LAR, AND TRAPEZOIDAL PRESSURE DISTRIBUTIONS 


Considering a laboratory sample, the equation for 
the exponent N is 


act \ 
» 3) 


When a definite numerical value is assigned to q,, N 





has a numerical value, designated’as Ni. The expres- 
sion then becomes 
_™G) 

SES AR NA CS Sein ery ee (16) 


Co Lee 


t=period of load application, in minutes. 
d,=thickness in centimeters that the labora- 
tory sample would have if it were possi- 
ble for the voids ratio to equal zero. 
d,j—=twice the maximum distance water travels 
when flowing from a sample of d thick- 
ness sandwiched between two filters. 
c=coefficient of consolidation of compressible 
soil, in centimeters per minute. 


To express ¢ of equation (16) in years instead of min- 
utes, and d) in feet instead of centimeters, the following 
computations are necessary: 


where 











1 centimeter= 5-55 feet 
and 
1 minute= : years 
365.24 «24 « 60 ; 
Then 
m(S) 2 
(2. 04X 1 2.54% 12)? | u(Z y 
DAT On Ab 24900460 a1 40002 


In this equation d, is in feet. 

For a compressible stratum in the field, sandwiched 
between two permeable layers of soil, the time required 
for a given percentage of consolidation to take place is 
computed from the equation, 


ni(3') 
ee sprang a ES a (17) 











where D,—the thickness in feet of the stratum at e=0. 

When a compressible soil stratum rests on rock and is 
covered with a porous layer, and when the pressure dis- 
tribution is uniform, the time-consolidation relation is 


When a soil stratum rests on impermeable rock and 
is covered with a porous layer, and the pressure dis- 
tribution is triangular (case 2), the expression becomes 


where N, is a definite numerical value corresponding 
to an assigned value for g:, the percentage consolida- 
tion at the time, ¢, for this type of pressure distribution 
and boundary conditions. For values of q: and q2, the 
corresponding numerical values of N; and .V; are given 
in table 2. 

For triangular pressure distribution (case 3) with a 
permeable layer only at the upper bounding surface of 
the compressible stratum, the time-settlement relation is 


For values of qg; the corresponding numerical values of 
N, are given in table 2. Here g3; denotes the percentage 
of consolidation at any time, t, with this type of pres- 
sure distribution (case 3) and with the stated boundary 
conditions. 

The procedures followed in computing the theoretical 
values of N,, N2, and N3, corresponding to values of 
91, 92, and q3 will be shown on pages 7 and 8 preceding 
the illustrative examples. 


TaBLE 2.—Values of Ni, Nz, and N3 corresponding to values of 
diy G2, and q3 for cases 1, 2, and 8, respectively 

















i 
Numerical values of— 
Consolida- 
tion q, 92, 
OF gs M Np Na 
Percent 

0 0. 00 0. 00 0. GO 

5 - 005 . 06 . 002 

10 . 02 4 005 
15 . 04 .18 01 
20 . 08 .25 02 
25 12 7 od 03 
30 pally nh) . 06 
35 . 24 .47 . 09 
40 aru . 55 12 
45 . 39 . 63 cally 
50 .49 73 . 24 
55 . 59 . 84 .32 
60 wile 95 . 42 
65 . 84 1, 08 . 54 
70 1. 00 1.24 . 69 
75 1.18 1, 42 . 88 
80 1. 40 1, 64 1. 08 
85 1.69 1.93 1. 36 
90 2. 09 2.35 Med 
95 2. 80 3. 02 2. 54 

100 Infinity Infinity Infinity 








With trapezoidal pressure distribution, the pressure 
resulting from the applied load may either increase 
(case 4) or decrease (case 5) with depth, and has a 
finite value at both the upper and lower boundaries or 
surfaces of the compressible stratum. The ratio of the 
pressure at the upper boundary to that at the lower is 
designated as vu. The value of u is unity for uniform 
and either zero or infinity for triangular pressure dis- 
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tribution. For trapezoidal pressure distribution, (case 
4), wis less than 1 and greater than zero. An interpo- 
lation factor, J, is applied in this case in computing the 
rate of settlement if the compressible stratum 1s overlaid 
with sand and rests on rock. When this condition 
obtains, interpolation is between the two extremes of 
uniform and triangular pressure distributions. With 
the compressible stratum bounded below by sand and 
above by rock, and with a trapezoidal pressure distribu- 
tion such that wu exceeds unity, (case 5) an interpolation 
factor, J’, is used in computing the rate of consolidation. 

The expressions for J and J’ in terms of wu are as 








follows: 
log| * eT 
For case 4, A= s ee Lek Ace Sag a (21) 
log— 
TT 
log] u(r—2)+2 
Horicasé-5, ~ J’= fue) UE Se PS 


Tv 


Values for J and J’ corresponding to values assigned to 
u are listed in table 3. 

The time-settlement relations for 
trapezoidal pressure distributions, are: 


cases 4 and 5, 


N.D,* 
Case 4, ‘=F 400¢ VOR See ees (23) 
N;D¢ 
Case 5, t= 400¢ yearsege seats (24) 
where PalVeew iN of GUNG SHV) ees. de J. ko (25) 
and N;=N;+ J" (N, —N3) SE SS ae (26) 


TABLE 3.—Values of u, J and J’ used in computing the rates of 
settlement for the two types of trapezoidal pressure distribution 




















Pressure increasing with Pressure diminishing with 
depth, uw varying from 0 depth, wu varying from 1 
to 1 to infinity 

Value of u Value of J Value of u Value of J’ 
0.0 1. 00 1.0 1. 00 
IN, . 92 
yf . 84 1.4 . 86 
1.6 . 80 
47) . 69 1] 1.8 oto 
1.9 03 
.3 OT 2510) ail 
2.2 . 67 
i . 46 225 . 62 
3.0 ait 
5 velo i) . 50 
210% 45 
6 ser 4.5 42 
5.0 39 
7 19 | 6.0 84 
7.0 . 30 
8 riage eA 8.0 27 
| | 9.0 25 
9 . 06 10.0 BY} 
| 12.0 . 20 
1.0 . 00 | 15. 0 17 
20. 0 13 
100. 0 02 




















FAIRLY CLOSE AGREEMENT FOUND BETWEEN THEORETICAL AND 
TEST VALUES 


The time-consolidation curve (fig. 2) represents the 
average curve for pressures of p=1.079, 2.118, 4.196, 
and 8.352 kilograms per square centimeter. The time- 
consolidation curve for the load p=0.61 kilogram per 
square centimeter was obviously influenced by experi- 
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mental error and therefore was not included in com- 
puting the values for the average curve. 

The divergence of the average curve from the curve 
given by equation (15) may be determined by means 
of equation (16) as follows: 

For different values of g,, say 25, 50, and 75 percent, 
the values of t are obtained from the curves. The 
average values of e at the different stages of loading 
are also obtained from these curves. Then the values 
of c, corresponding to the average values of e, are read 
from the curve in figure 3; the value of d) is computed; 
the numerical values are substituted in equation (16); 
and the values for N, may thus be computed, assuming 
that the relationship in equation (16) obtains. 

The thickness of the sample, d, at a voids ratio, 
e=1.670, was found in the test to be 1.130 centimeters. 
Therefore the thickness d) when e=0 would be 0.423 
centimeter. 

The values for N, thus obtained (see table 4) for 
equal values of q are averaged and found to be 0.11 
for q:=25, 0.44 for q=50, and 1.30 for q,=75. 

The theoretical values of N, for the same values of 
g, are shown in the last column of table 4. The agree- 
ment between the determined and theoretical values is 
fairly good. The procedure in computing these theo- 
retical values of N, follows: 


TABLE 4.— Values for qi, e, c, and N, obtained from laboratory data 


























- .¢ | Walues of 
ae Values of c | Values of 7 
Pressure Consoli- ae ae correspond-| Ni com- | ~ Lee 
interval dation, | patio. ¢ ing to puted om 
n AES) average € | from data!) theory 
Kg per sq.cm | Percent 
{ 25 lea 0. 0108 0.10 0.12 
0.610-1.079__-_ 50 1.138 . 0116 .49 .49 
| rH 1. 08 . 0140 1.47 1.18 
| 25 1. GO . 0152 18 12 
1.079-2.118__-_ 50 . 97 . 0159 .42 .49 
| TB . 94 . 0166 1. 20 1.18 
25 . 87 0177 me a1. 
2. 118-4. 196___ 50 . 84 . 0180 . 42 .49 
h3) 80 . 0183 1,32 {ats 
| 25 . 74 . C188 . 09 w12 
4, 196-8. 352___ 50 72 . 0189 . 44 49 
| 75 69 . C191 1. 22 1.18 
| 











1 For q:=25 percent, average Ni (computed from data in column 5)=0.11. For 
qi=50 percent, average Ni (computed from data in column 5)=0.44. For qi=75 per 
cent, average Ni (computed from data in column 5)=1.30. 


COMPUTATIONS GIVEN FOR VARIOUS TYPES OF PRESSURE DIS- 
TRIBUTION AND BOUNDARY CONDITIONS 


Case 1.—The computations for case 1 cover uniform 
pressure distribution, irrespective of boundary condi- 
tions, and triangular and trapezoidal pressure distri- 
bution when the compressible stratum is bounded on 
both top and bottom by permeable strata. To compute 
the value of N, corresponding to any definite value as- 
signed to g; in the expression, 

n=1—S(o¥ + ge +5504 Se dane .({15) 


1 


one may, as an approximation, neglect all terms of e 


except the first. Then 
ate or 
Sl Saxer sy ee ge ee (2 
qg=1 ue -------- (27) 
or 
8 f (JQ 
1—gi=—40°=--------------- (28) 
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For q=s5» equation (28) becomes 


Oe 
pe 1 
and 
1 8 7 
log e 5 log e pos: 


) 


where “log .”’ designates the natural logarithm. 


SP ee 050 lomo lone 
N,=log « 5—log , 5=2-8026( log ~—log 5) 


where “‘log’’ denotes the base 10 system of logarithms. 
Hence NV,=0.48. 

The more accurate procedure would be to assign 
values to N, include several terms of the series, and 
solve for g;. In this manner it is found that for N,;=0.49 
the value of gq, is 0.50 or 50 percent. (See table 2.) 
However, it is seen that by neglecting all terms of e 
except the first in equation (15), the error, 0.49—0.48 
=0.01, is small. 

Case 2.—Triangular pressure distribution with zero 
pressure at the upper boundary, when the compressible 
stratum is bounded by an impermeable stratum below 
and a permeable stratum above. For this condition 
the theoretical time-settlement relationship is given by 
the expression, 

32 1 


Q@=1-S (e-N— Aen tytn sie a 


re 27 aR ATos Tee) 


By neglecting all terms of e except the first, one may 
write, 


32 
1, acy * Pe es By ee St | (30) 
; 1 ‘ 39 
for =~ 5’ N »=log=3—log.5=0.72. 


This compares favorably with the value, N.=0.73 for 
q2—50 percent, table 2, obtained by taking more terms 
in the series than the first in equation (29). 

Case 35.—Triangular pressure distribution with zero 
pressure at the lower boundary, when the compressible 
stratum is bounded by an impermeable stratum below 
and a permeable stratum above. It will be shown in 
part 2 of this paper that 


d3=291—e a a Ni ee -(31)} 


The simplest procedure in computing values of Ny; that 
correspond to values of gs is to take the case, N;=N, 
=N;3. For Ni=N,=N3=0.12, it is seen in table 2 that 
g=0.25 and @=0.10, so that g=2<0.25—0.10=0.40, 
That is for g,=40 percent, N;=0.12. 

Jase 4.—Trapezoidal pressure distribution with pres- 
sure increasing with depth, when the compressible 
stratum is bounded by an impermeable stratum below 
and a permeable stratum above. The value, Nj, is 
computed from equation (25), which is Ny=N,+ 
J (N.—N,), applicable for m=qa=q. J in turn is 
computed from equation (21) for any given value of w. 
For example, when uw=0.50, J=0.36 (see table 3). For 
1=d@=q=50 percent, N,=0.49 and N.=0.73. Then 


N,=0.49+0.36(0.73—0.49) =0.58. 











When this value is substituted for N, in equation (23), 


_ ND? 
~ 1,400 ¢ 





t years, 


the time required to produce 50 percent of the ultimate 
consolidation is easily calculated. The value of ¢ is 
obtained from laboratory data according to the proced- 


ure already indicated and D)= as already shown 


dD, 
1 + Al 
(equation 2). The thickness, D,, of the compressible 
stratum prior to applying a load at the surface is known 


from the test borings. If D,=20 feet and e,=1.5, then 


Tee feet. And if c=0.01 and N.4=0.58, then 


Aes 


0.58 X 64 
14 


i= This is the time required for 


50 percent of the ultimate consolidation to take place 
under the given conditions, since a value of N4y=0.58 
corresponds to 50 percent consolidation. 

Case 5.—Trapezoidal pressure distribution with pres- 
sure decreasing with depth, when the compressible 
stratum is bounded by an impermeable stratum below 
and a permeable stratum above. For q=q=q;, the 
value, N; is computed from equation (26), which is, 


N;=N3+d’ (Ni—Ns3), 


for any given value of wu. From table 3, J’=0.71 when 
u=2. For 1=¢=¢;,=—50 percent, N;=0.49 and N3= 
0.24 from table 2. Then by substitution, N;=0.24+ 
0.71(0.49—0.24)=—0.42. Now take the conditions al- 
ready assumed for case 4, that is, c=0.01, e,=1.5, the 
average voids ratio in the stratum prior to loading, 





=2.65 years. 





and D=20 feet. Then 
t= ND sp Uae r 1.9 years 
TT AQ0E) tee ee 1 teat he ee 


DISCUSSION OF CONDITIONS GOVERNING ILLUSTRATIVE EXAMPLES 
THAT FOLLOW 


In the following examples it is assumed that the 
material composing the compressible stratum has the 
consolidation characteristics derived from computations 
of the data plotted in figures 2 and 38, and that the 
weight per cubic foot of this material as well as that 
of the sand above it is 100 pounds. Cases 1 to 5 
inclusive are linear types of pressure distribution, that 
is, the pressure distribution diagrams are composed of 
straight lines. All of these five types of pressure dis- 
tribution give the same time-consolidation relation, 
equation (17), when the compressible stratum is 
bounded both above and below by drainage courses 
(boundary condition A). Equation (18) is the time- 
consolidation relation when there is but one drainage 
course (boundary condition B) and the pressure distri- 
bution is rectangular (case 1). When there is but one 
drainage course, the stratum resting on rock and being 
overlaid with sand, equations (19), (20), (23), and (24) 
are the time-consolidation relations for cases 2, 3, 4, 
and 5, respectively, and computations are made with 
the aid of tables 2 and 3 and the interpolation formulas, 
making use of J and J’, that is, equations (21), (22), 
(25), and (26). 

Cases 1 to 5 inclusive may be defined in terms of w, 
the ratio of the pressure at the upper surface of the 
compressible stratum to the pressure at the lower sur- 
face, that is: 
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Case 1, rectangular pressure distribution, uw—1. 

Case 2, triangular pressure distribution, w=0. 

Case 3, triangular pressure distribution, w—infinity. 

Case 4, trapezoidal pressure distribution, w varies 
from 0 to 1. 

Case 5, trapezoidal pressure distribution, w varies 
from 1 to infinity. 

In solving any problem involving case 4, the formulas 
and principles of cases 1 and 2 are applied. When the 
pressure distribution is of the type case 5, application 
is made of the formulas and principles that pertain to 
cases 1 and 3. Cases 1 and 2 are both limiting condi- 
tions of case 4, and cases 1 and 3 are both limiting 
conditions of case 5. The types of pressure distribu- 
tion most frequently encountered in foundation prob- 
lems are cases 1 and 5, more especially 5. 

Example 1.—Case 1 pressure distribution; type A 
loading. Consider a uniform load of 3 tons per square 
foot over an area of great extent. Immediately be- 
neath the loaded area there is a layer of sand, 10 feet 
thick. Below the sand is a compressible stratum also 
10 feet thick and below this stratum there may be 
_ either (A) sand or (B) rock. It is desired to know the 
total settlement and the time required for 25, 50, and 
75 percent of the total settlement to take place. 

Use equation (17) for boundary condition A and 
equation (18) for boundary condition B for the time- 
consolidation part of the problem. It is assumed that 
consolidation of the compressible stratum resulting 
from its own weight and that of the sand above it has 
reached completion prior to application of the surface 
load. 

Prior to loading, the pressure is 1,000 pounds per 
square foot at the top and 2,000 pounds per square foot 
at the bottom of the compressible layer. In figure 2, 
these pressures correspond to voids ratios of 1.28 and 
1.06 respectively. A pressure of 1 kilogram per square 
centimeter is very nearly 1 ton per square foot so that 
the units of pressure of figure 2 (load-compression curve) 
are directly applicable. The average voids ratio prior 
to loading is then 


1.28+1.06 _ 


”») 


“a 


alee 


After loading, the pressure at the top of the stratum is 
3.5 tons and at the bottom, 4 tons per square foot. 
These pressures correspond to the voids ratios (see fig. 
2), 0.80 and 0.78, respectively, the average, e, being 
0.79. The total settlement is computed by equation 
(14), that is, 


65 1.17—0.79 
J be, 1+1.17 


If there were no voids whatsoever in the material 
composing the compressible stratum, its thickness 
would be 


Dy 


Sb X10=1.75 feet. 





1 1 
mea 30.17 


The average voids ratio is reduced by the surface load- 
ing from ¢;=1.17 to ¢=0.79. Within this interval 
‘ ed fas 
there is an average ¢ of Sete elo: 
a voids ratio of 0.98 corresponds to a value for c, the 

coefficient of consolidation, equal to 0.0157. 





<10=4.61 feet. 


From figure 3, 











Now apply equation (17) for boundary condition A, 
the compressible stratum sandwiched between two 
drainage courses. The times required for 25, 50, and 
75 percent of ultimate consolidation (1.75 feet) to take 
place are, using table 1 for values of N,, 


(2) 


for q,=25 percent, and V,=0.12 AES a 
Ssh a p and NV, =0.12, 1,400 ¢ 


0. i2x(4*) 


1,400 X 0.0157 





=(.03 year or 11 days, 
0.49 43") 
* @,=50 percent, N,=0.49, t=— + = 4 
for ¢:=50 percent, N,=0.49, ¢ 1,400 0.0157 


6 z 
Lasx(8) 


=0.12 year=44 days, 


*? 
1,400 X0.0157 
=(0.285 year=104 days. 


p 


{Ol Gee (> percent. avi == 1.18, -0=— 





For condition B, the stratum covered with sand and 
resting on rock, the values of N, corresponding to 
values of g; are the same as for condition A. Solving 
0.12 (4.61)? 9 45 
400001 bv ae 
year=44 days. Similarly, for q:=50 percent, t=0.48 
year, and for q,=75 percent, t=1.14 years. 

Twenty-five percent of 1.75 feet, the ultimate settle- 
ment, is 5.2 inches. For condition A, this amount of 
settlement occurs in 11 days; for condition B, 44 days 
are required for this amount of settlement. 

Example 2,—Case 5 pressure distribution; type C 
loading. Instead of a loading area of great extent as 
in example 1, there is a rectangular rigid footing, 10 
feet wide and 30 feet long, carrying a uniform load of 
3 tons per square foot. This is at the ground surface 
above the same compressible stratum described in 
example 1. The average voids ratio, e,, of the com- 
pressible layer prior to loading is taken as 1.17; its 
thickness prior to loading is 10 feet; and the weight of 
the sand at the upper level of the stratum is 1,000 
pounds per square foot. 

The vertical distance, z, of the footing above the 
upper boundary of the stratum is 10 feet. The ratio 
of the length of the footing to zis 3. This ratio is too 
large to warrant taking the uniform load on the footing 
as a point load at its center. For purposes of computa- 
tions it 1s necessary to consider the footing as sub- 
divided into rectangles such that the ratio of the length 
of a single rectangle to z does not exceed one-half (see 
figs. 5 and 7). The footing may be divided into twelve 
5-foot squares to meet this requirement. 

The total load on the footing is 310 30=900 tons. 
The load on a single square is 75 tons and it is con- 
sidered that this load is concentrated at the center of 
the square. The horizontal distance from the center 
of a square to the center of the 10- by 30-foot{foundation 
is r and this is found from the properties of a right 
triangle. For example, the center of square number 5 
from C, the center of the footing, is the square root of 


yi tar 2 
@) +(3) sor Bo" The same value for 7 is obtained 





equation (18) for q,=25 percent, t= 
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Figure 7.—REcTANGULAR LOADED AREA (TyYpb C). 
7 — see ee 
for squares 6, 7, and 8. The ratio, — appearing in 

& 


: ; Bs) l 5 

equation (10) is ———= or ~—= for squares 5, 6, 7, and 8. 
10/2 242 

These values and those obtained for the other squares 

(fig. 7) are given in table 5, together with the value of 


<td ; ate ihe 

kK obtained by substituting the value for — in equa- 
& 

tion (10), 


K=— xy 


The three values for A, corresponding to the three 
groups of load numbers given in the first column of 
table 5, are added together, their sum being 0.5404. 
This must be multiplied by 4 to give K for the entire 
footing. The vertical stress at a point 10 feet below 
C and on the vertical center axis through C is then 
calculated from equation (9). That is, 

75 


?.=—> 100 <4 «0.5404 = 1.62 tons per square foot. 


TaBLE 5.—Computations of vertical stresses below a rectangular 
footing 10 by 80 feet 
































Load numbers (see fig. 7) Distance, r Ratio, = Been 
Feet | 
5 I 
5, 6, 7, and 8- > CO Oba =e. 354 0. 3555 
| ees A Ree | 
| 3, 4,9, and 10. . ‘ bended - = 0.790 . 1420 
: a 5v26) |) 26. : 
1, 2, 11, and 12 | 9 | + =1.274 0429 
Total | 5404 
By taking z=20 feet and computing in the same 


manner, it is found that at the lower boundary of the 
compressible stratum and on the vertical center axis of 
the footing, p, is 0.73 ton per square foot. The pres- 
sure is & maximum for points along the vertical center 
axis of the footing and the distribution of this maximum 
pressure is trapezoidal, case 5. At points removed 
from the vertical center axis there is also case 5 pres- 





10 BY 30-FT. RECTANGULAR FOOTING 


Sc ESaND TS Gaee 

























3,240 


ACTIVE 





1460 


ie CANO tee 


Li 
== 7, )000,10007 o> 'yan yas 


Ficure 8.—DistriputTion or Maximum PRESSURE UNDER 
REcTANGULAR Footine. Sort LAYER CONTAINED BETWEEN 
Two PERMEABLE SAND LAYERS. 


sure distribution, but the value of p, at any level is 
less than the value of p, for the same level on the 
vertical center axis. For computations consider only 
the maximum pressures, which are along the vertical 
center axis. 

The shaded trapezoidal area marked ‘inactive’ in 
figure 8, denotes the pressure distribution resulting from 
the weight of earth material alone. The trapezoidal 
figure marked ‘active’ denotes the distribution of max- 
imum pressure in the compressible stratum that results 
from the uniform load on the footing. It is again 
assumed that consolidation of the stratum caused by 
its own weight and that of the sand is complete at the 
time the surface load is applied and hence all of this 
pressure has been transferred from water to soil and 
therefore cannot act in combination with the applied 
surface load in causing consolidation. The applied load 
producing the active trapezoidal pressure distribution 
produces further consolidation beyond that already 
effected by the weight of earth material. 

When the surface load is applied, the pressure per 
square foot at the upper boundary of the compressible 
layer is 1,000+3,240=4,240 pounds or 2.12 tons. That 
at the lower boundary is 3,460 or 1.73 tons. From 
figure 2, these pressures correspond to voids ratios of 
0.90 and 0.94, respectively, the average value of e being 
0.92=e. Prior to loading, the average value of e, as 
found in example 1, is 1.17. The ultimate settlement 
is therefore 

722 ().09 
= “5 aX 10=1.15 feet=13.8 inches. 


The average voids ratio is reduced by the surface loading 
from 1.17 to 0.92. The average of these is 1.045 and the 
corresponding value of ¢ from figure 3 is 0.0138. 

Example 2 (continued).— Boundary condition A; the 
stratum is sandwiched between sand layers. For the 
time-consolidation computations, again apply equation 
(17) for boundary condition A (see fig. 8). 


yy (Poy A GIN 
EN OS eg 12 Ne 


1,400 ¢ 1,400 0.0138 














For g,=25 percent, t= =0.033 


year=12 days. 
Similarly, =0.135 year=49 days for ¢,=50 percent and 
0.325 year=119 days for g,=75 percent. 
These values are nearly equal to those computed for 
boundary condition A, example 1, and would be the 
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same were it not for small variations in ¢ as the pres- 
sure is varied. For most practical purposes, ae 
boundary condition A obtains, the time-consolidation 
relation can be considered independent of the amount of 
pressure and the type of pressure distribution. 

Example 2 (continued).— Boundary condition B; the 
compressible stratum rests on rock and is overlaid with 
sand. When condition B obtains in this example (see 
fig. 9), use equations (24) and (26), together with tables 
2 and 3. The total settlement, Q, is the same for both 
boundary conditions, A and B. The value for uw (con- 
sider only the active pressure, fig. 9) is 2.22. Then 
J’ is obtained from table 3 and is taken as 0.67. From 
table 2, when q=q3=q;=25 percent, N,=0.12 and 
N;=0.03, and from equation (26), 


N.—N,-+.J’ (N,— N;)=0.03 +0.67(0.12—0.03) =0.09. 


This value for N; 1s substituted in equation (24) and 
taking c=0.0138, then for g;=25 percent, 








NED ADS ICS A) 
PemOte 00 001ss oy ot 
For 4=4¢3:=¢q;=50 percent, N, (table 2)=0.49 and 

INo== 0.24. 


Then N;=0.24+0.67 (0.49—0.24)=0.41 and 


0.41 (4.61)? 
1,400 0.0138 


Similarly for g;=75 percent, 


j= 1.08 X (4.61)?_ 
1,400 < 0.0138 


[p= 





=0(0.451 year=165 days. 


=1.19 ye 435 days. 





Example 3.—Case 4, pressure distribution; type E 
loading; boundary condition B. It would seem that 
case 4 pressure distribution is a comparatively rare 
occurrence. Problems of this type may be encountered, 
however, when consolidation produced by the weight of 
earth material alone is inc omplete at the time of load 
application at the ground surface. Figure 10 illustrates 
such a condition. If it is assumed that consolidation in 
the compressible stratum 5 feet thick resulting from the 
weight of sand above is only 50 percent completed, then 
this weight exerts an active pressure which is zero at the 
top and ai ,000 pounds per square foot at the bottom of the 
compr essible layer. The distribution of active pressure 
resulting from the weight of sand alone is, therefore, 
triangular ; 





It is ened that consolidation of the 
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ReEstTiInc ON Rock. 


compressible layer as produced by its own weight is 
complete at the time of loading. 

Referring to figure 10, a load of 3 tons per square foot 
is applied over a rigid circular footing 10 feet in diameter. 
It is desired to determine the total settlement and the 
rate of settlement of the compressible layer which rests 
on rock and is covered with the sand layer 10 feet thick, 
50 percent of the consolidation producible by the weight 
of sand being completed at the time of load application. 

To compute the pressure, p,, at a point on the vertical 
centerline 10 feet below the rigid footing, substitute 
values in equation (12). The cosine of angle 8 (fig. 10) 








a LO 2 
is ~—=—= 7=and p equals 3 tons per square foot. Then 
i195 
V125. 5 
5 
P2=p (1—cos®Bs) =; (1-55;)=0. 854 ton per square 
Bale 
foot. 


At a point 15 feet below the footing and on the verti- 
cal center axis, p,=0.438 ton per square foot. There- 
fore, the applied load produces a pressure of 1,707 
pounds per square foot at the upper boundary and 876 
pounds per square foot at the lower boundary of the 
compressible layer. 

Prior to loading, e at the upper surface of the layer is 
1.28, corresponding to a load of 1,000 pounds per square 
foot, and at the lower boundary it is 1.44, corresponding 
to a load of 500 pounds per square foot. The average 
is 1.36=e,. The voids ratios (see fig. 2) corresponding 
to pressures of 1,000-+-1,707 pounds per square foot at 
the top and 500+1,876 pounds per square foot at the 
bottom of the layer are found (fig. 2) to be 0.99 and 1.01, 
respectively. The average is 1.00=e. The total set- 
tlement is, therefore, 


ou ie 36—1.00 
1+1. eon 

It is important to note that prior to loading, soil at the 
top of the compressible stratum has reached equilibrium 
under the weight of sand (1,000 pounds per square foot) 
and that soil at the bottom of the stratum has not as yet 
been consolidated by the sand but has been consolidated 
by the weight of compressible soil, 1. e., 500 pounds per 
square foot. Therefore, the voids ratio at the top of the 
layer prior to loading corresponds to a pressure of 1,000 
and that at the bottom to a pressure of 500 pounds per 
square foot. After a very long period of time subse- 
quent to loading, the soil at the top of the compressible 
layer will be in equilibrium with a pressure of 2,707 
pounds per square foot and that at the bottom with a 


5=0.76 foot. 
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STRIP OF INFINITE LENGTH are the same as in example 3 except that in this case 

(OFEET WIDE the rigid footing is a loaded strip 10 feet wide and of 

Ae relatively great length instead of a circular bearing area 

PER SQ.FT, | and consolidation caused by weight of sand is complete 

BEER z Say AG prior to loading. Figure 11 shows a vertical cross 
i Bees. ae i section. To find p, at the top of the compressible 
Bote tera oa - ogee Gree layer and at a point on the vertical center axis through 
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Sor, LAYER OVERLAID 


pressure of 2,376 pounds per square foot. These pres- 
sures correspond to definite voids ratios, found experi- 
mentally (see fig. 2), and are used in computing Q. 
When load is applied, the active pressure at the top 
of the layer is 1,707 and that at the bottom is 1,876 
W707 
miko a 
value between 0 and 1 (case 4). Apply equations (23) 
and (25) and use the values for V, and N, corresponding 
to different values for q;=q, table 2. Values of J are 
found from table 4, and for w=0.91, J=0.06. For 
=G@=25 percent, Ni=0.12, N.=0.31, and 


pounds per square foot. Hence wu 


N,=N,+d(N2—M,) =0.12-++0.06(0.31—0.12) =0.13 





1 1 
ees -X5=- =, X5=2.12 fe 

Dy 1 pene 7.36 ) wil feet. 

For e=“ E “21.18, e=0.0106 (fig. 3). 





4 : IND) | OMS A212) 

a . —P95 A == = eos = 
PO De TOUE Gi 25 (DELCeNts te ernne man O UI On 
0.039 year=14 days. 

0.50 (2.12)? 








Reet acces | emule =F ‘cent == = fs 
Similarly, for q;=50 percent, ¢ 1,400 <0.0106 0.15 
year=55 days, 
EES COPA 
and for qs=75 percent, t= aes ) =(0.36 year= 


1,400 0.0106 
132 days. In 132 days the settlement is 75 percent 
of 0.76 foot or 6.8 inches. 

Example 4.—Case 5 pressure distribution; type D 
loading; boundary condition B. All other conditions 





the cross section, apply equation (11). The angle, a, 
is expressed in radians and sin a=0.80. 
0,097 (0.93 +0.80) =3,290 pounds per 


p.—"(a-+sina) = 
square foot. 
At the lower boundary 
p= 7 (0.64-+0.60) =2,375 pounds per square foot. 
This is case 5 pressure distribution, and w=1.38. The 
average voids ratio prior to loading is the average of ¢ 
corresponding to 1,000 and of e corresponding to 1,500 
pounds per square foot and is 1.21=e,. Subsequent to 
loading the values of e at the top and bottom of the 
compressible stratum correspond to the equilibrium 
pressures, 4,290 and 3,875 pounds per square foot, 
respectively. From figure 2 these values of e are 0.90 
and 0.92, the average being 0.91=e,. The total settle- 
ment is therefore 


and on the centerline, 





.21—0.9 i 
Q-|+ 5 = * [x5=0.68 foot or 8.2 inches. 


To compute the rate of consolidation, apply equations 


(24) and (26) and make use of tables 2 and 3. J’=0.87 

| for uw=1.38 (by interpolation). For peed 55 = 
0.0134. 

N;Di 0-126 ae 








= = —% j TC = =F gar 
HOD Uso dem: 2) VeTCen inc eee conan 


0.03 year=11 days. 
0.46 X (2.26)? 
1,400 0.0134 





For g;=50 percent, {= =(0.125 year=46 


days. 
1.14 X (2.26)? 
1,400 X 0.0134 





For ¢;=75 percent, t= =0.310 year=113 


days. 
In 113 days the strip settles 75 percent of 8.2 inches 
or 6.2 inches. 

The preceding examples are intended to illustrate the 
application of the theory of soil consolidation. Part 2 
contains the derivation of equations and the develop- 
ment of the basic theory. 
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PART 2.—DEVELOPMENT OF THEORY OF SOIL CONSOLIDATION 


In 1885, J. Boussinesq ® developed in considerable 
detail and usable form a solution of the general problem 
of determining the distribution of stresses and strains 
in a semi-infinite, elastic, isotropic solid, bounded by a 
plane and loaded by a single concentrated force at a 
point on that plane. The formulas of Boussinesq have 
had wide application in the development of the science 
of soil mechanics. 

The distribution of stresses across any plane parallel 
to the plane surface is shown in figure 12. If a single 
concentrated vertical load, P, is applied at a point on 
the surface, then at a radial distance, R, from this 
point and at a depth, z, below the surface, the stresses 
at the point (x, y) are: 


SLaNe ; - 
Pr= 3 RP” parallel to R 2 (32) 
Hele hrd : 
Pio 5 7 pe fae parallel to x axis wine) 
areey- 
Pu 37 fe parallel to y axis. _ (34) 
anes 
P= 57 Pe Pe parallel to z axis_ Sb) 


The derivation of these equations may be found in 
almost any standard text dealing with the theory of 
elasticity.’ These equations contain no elastic con- 
stants, however, and for this reason it has probably 
been assumed that the distribution of stresses does not 
depend upon the type of material. Cummings,® 
Krynine,’ and others have discussed this point at length. 

A complete analysis of the stress distribution includ- 
ing planes not parallel to the boundary involves 6 com- 
ponents of stress, 3 shears and 3 normal stresses, and 
the elastic constants appear in some of these equations. 
In the theory of consolidation of soils we are concerned 
mainly with the vertical normal stress, p,. At a given 
depth, z, it is seen from equation (35) that p, is greatest 
when F is smallest, that is, the maximum values for p, 
are along the vertical centerline of the loaded area 
e—y—0; fie. 12). 

The value of p, at any point is easily computed, 
knowing that R’?’=r?+z2?=27?+y?+27. If a combi- 
nation of point loads is given, we may obtain the stresses 
at a point resulting from each load considered sepa- 
rately, then add these and obtain the total pressure at 
the point produced by the combination of loads. 
Gilboy 1° has shown that by dividing a uniformly 
loaded area into rectangular elements so that the 
longer side of any element is less than half the distance 
from the center of the area to the point, the computed 
value for p, at the point will not deviate from the 
correct value by more than 6 percent. The accuracy 
is greater with further subdivision of the rectangular 
elements. 


6 Application des Potentiels a l/Etude de l’Equilibre et du Movement des Solides 
Elastiques. Paris, 1885, ee ; 

7 For example, see chapter 11 of Theory of Elasticity, by S. Timoshenko. F 

8 Distribution of Stresses under a Foundation, by A. E. Cummings, Proceedings, 
American Society of Civil Engineers, vol. 61, no. 6, August 1935. ' 

* Discussion of Cummings’ paper by Clement C. Williams, D. P. Krynine, and 
L. C. Wilcoxen, Proceedings, American Society of Civil Engineers, vol. 68, no. 8, 
pt. I, October 1935. : : ; 

10 Barths and Foundations. Progress Report of Special Committee, Proceedings, 
American Society of Civil Engineers, vol. 59, no. 5, May 1933. 
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Figure 12.—DiIsTRIBUTION OF STRESSES ACROSS ANY PLANE 
PARALLEL TO THE PLANE SURFACE. 


FORMULAS DEVELOPED FOR DETERMINING PRESSURE AT POINT 
BELOW LOAD 


With reference to figure 12, if the expression for 
a) 
Gea eS is integrated in the y direction from 
=—o to y=+o and the resulting ‘line load”’ is 
multiplied by dz, we obtain the total load on a strip 
infinitely long and of width, dr, produced by the point 
load, P, at the boundary. In this integration, x and z 
are constant. Thus: 


: LAER 3P2 
l 2 QrR> Q2a(a?4+ y?+ 2?)5? 


3 P23 (°° x? 21 p27-6/2 n Pe ec 





This is easily integrated if we put =tan 60 


Y 
7? 4 ep ak! 
(see fig. 12). Then 


83P2 6=7/2 : 
\z= (dr) rs | cos*6d6, 





or 


ore ; sin’9 ]*/? IPS 
Jo saegapy sin 0 [tem aap de 80) 


If this expression is again integrated in the x direc- 
tion, taking twice the value obtained between the limits, 
z=0 and c= © (fig. 12), the result is equal to P, the 
applied force. Obviously, on any plane of infinite 
extent and at a finite distance below and parallel to 
the loading surface, the total distributed load is equal 
to the applied load at the surface. The mathematical 
result js in agreement with this self-evident fact. 
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LoapEp AREA oF INFINITE LENGTH AND OF WIDTH 2b. 


In equation (36), if P is taken as p, the ver tical pres- 
sure per unit area, and if p is considered as concentrated 
at the center of a unit area, then the total vertical 
pressure on. Re strip at a distance below the surface is 





equal to ye de. This is the total pressure, /p., 


T(x mare 
on a strip below the surface produced by a unit load on 
the surface. Conversely, the total load on an infinitely 
long strip at the surface and parallel to the y axis, 
produced by applying the unit load, p, over the length 


23 
of this strip will cause a vertical pressure —3~— ae eee Sa 2 de r, at 
a point below the surface. Briefly, a point ise at the 


2pe 
ee 
below the surface and conversely, a strip load of unit 
load per unit area on the surface produces a pressure, 

2p2 
Ree) 
reaction being equal. It will be noted that the integra- 
tions of the Boussinesq formulas are with reference to 
areas or points below the loaded area at the surface. 
At the surface, z=0 and there is nothing to integrate 

that has reference to the fundamental equation (35). 
Figure 13 1s a vertical cross section of a uniformly 
loaded area (p per unit of area) of infinite length and 
of width, 2b. The load on a strip of infinitesimal width, 
dx, produces a& pressure, p,, at a point T equal to 


aDe2 
pee 
a(2?+ 2? 
limits, hee is obtained the value for p, at T resulting 
from the total load over the whole area, infinite in 


length and of width, 2 The integration is more easily 


“cot 6 (fig. 13). Then at the 


surface produces a total load, - aCe dx on a strip 


5dr at a point below the surface, action and 


If this expression is integrated between 


performed if we let 


S rTlN 
point T, 


2) me a 
De eal 1+- =). dex 72 | ( 


eee ee (37) 


1+ cot?6)~*esc76d6 


="[a—sin a COS(a,+ a2) ]_- 


If the point, T, is on the vertical centerline through 
the vertical cross section (fig. 13), then a2=180°—a, 
and equation (37) reduces to 


P:= E (qa: sina) - 
T 


for uniform pressure distribution (p per unit area) 
over a circular bearing area of radius r (see fig. 14) the 











AREA. 


Ficure 14.—Circutar LOADED 


pressure at a point below the loaded area and on the 
vertical centerline is, 


mee vf 3p ae ye 53 iE. : 
eel) lh De EL SPF |, Gaye 


Integrating and simplifying, 


(145) 


MECHANICS OF SOIL CONSOLIDATION OUTLINED 


and since — J=tang 





0p | 1 
Cz 


Consider the flow or egress of water from a wet, 
compressible, soil stratum subjected to vertical loading. 
The soil is relativ ely impervious and is “‘sandwiched’’ 
between two permeable sand-gravel layers, L, L, figure 
15. The movement of water saturating the stratum 
is in both the upward and downward directions, 
parallel to the vertical axis, ZZ, and directed toward 
the two permeable layers, L, L. 

At any finite time, t, after load application, that part 
of the total pressure ‘carried by the water, p,, 1s consid- 
ered as the “hydrodynamic excess.” " The remainder 
of the load is carried by the solid material, the ‘“‘soil 
skeleton.” At the outset of loading (t=0) all of the 
pressure is carried by water if the compressible layer 
is saturated. At any subsequent time, the hydrody- 
namic excess, Py», 1s zero at the surfaces of contact of 
the soil str atum and the permeable layers, but elsewhere 
it is a finite quantity, varying with the depth, z. The 
upper and lower faces, AB and CD of the elementary 
prism, ABCD, are each of unit area. The thickness 
of the elementary prism is Ag (fig. 15). 

We may then define the following expressions: 

p=the total pressure per unit cross-sectional area. 
It includes the pressure on the soil skeleton and the 
pressure, P,, carried by water alone. 

Pw=the pressure per unit cross-sectional area that is 
carried by water alone. 

P— Pw the pressure per unit cross-sectional area that 
is carried only by the soil particles, i. e., by the soil 
skeleton. 

a=the rate of change of the voids ratio with respect 
to (p—p») and is called the coefficient of compressibility. 

eae reduced coefficient of permeability; see equa- 
tion (4), Part 1. 


Ap» 
= “=the rate of change of the hydraulic pressure, 





Pw, With respect to the depth, z. 
pressure gradient. 


UD w 
Spee ==the rate of change of ae with respect to 2 and 


It is the hydraulic 





1 Terzaghi, K., Erdbaumechanik. Franz Deuticke, Vienna, 1925. 
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is taken as constant through the very small vertical 
distance, Az. 








LP Waihies SE one 
Az = =the total change of gE within the inter- 
val, Az. 
a? Pw : : , 
ky > 4ze=the rate at which the elementary prism 


loses ator It represents the difference between the 
quantity of water flowing from the upper face, AB, in 
unit time and the quantity flowing into the prism 
through the face, CD, in unit time. It is therefore 
the rate of crit reduction of the elementary prism. 

















dp d(p—p : 
i= (7 a tthe rate of change of p—p, with 
respect to time. 
dp , ar 
a-Fe= the rate of change in voids ratio with respect 
to time. 
abe ie 
ane z=the rate of change in volume of the ele- 
mentary prism with respect to time. 
dp» LD» dp, 
The expressions an = > and - - % are considered as 


the partial derivatives. The faces, AB and CD, of the 
elementary prism are both of unit horizontal cross- 
sectional area and are parallel to each other. 

The amount of water lost by the elementary prism 
LD 
dz 2 
tween the amount of water flowing out from the face 
AB and thet flowing into the prism through the face 
CD. But the rate of loss of water (which i is the rate 
of reduction in volume) by the elementary prism is also 
Za _ ihe 


in unit time is k,>~Az, which is the difference be- 








equal to aPers Then ko Bene ama Ae. Divide 
both sides of this expression by a and write ¢ for - os 
Then, since Az cancels, 

d? Pr _ _ Dw 

Cc dt 777777 (40) 


This is analogous to the Pest of the equation 
for the flow of heat. A more complete derivation of 
equation (40) has been published.” 


TIME-CONSOLIDATION RELATIONSHIPS FOR UNIFORM PRESSURE 
DISTRIBUTION DERIVED 
A solution of equation (40) is 


Dy»=e-~K*(A cos Kz+B sin Ka2)_. - (41) 


where A, B, and K are constants. 


12 Mathematical Theory of the Process of Consolidation of Mud Deposits, 
Alberto Ortenblad. Journal of Mathematics and Physics, vol. 9, no. 2. 


by 
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Figurm 16.—PRESSURE DISTRIBUTION IN A COMPRESSIBLE 
LAYER BETWEEN Two PERMEABLE LAYERS. AT ANY LEVEL, 
Z, Pw=THE PRESSURE TAKEN BY WATER. p=THE TOTAL 


PRESSURE. P—DP »=THE PRESSURE ON THE SOIL SKELETON. 
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Figure 17.—A, CONDITIONS IN COMPRESSIBLE 


Layer. At THE Instant Loap 1s APPLIED (t(=0), pw=p AT 
Aut Levers. B, PREssuRE CONDITIONS IN COMPRESSIBLE 
LAYER SUBSEQUENT TO LoAp ApPrLICcATION. WHEN ¢t>0, 


Dp OVA ee — (RANI ANe —= 20 ne yy <A AND LE ey BES 

In figure 16, the reduced thickness (the thickness the 
compressible stratum would have were there no voids 
between the particles) of the soil layer is 2/) where J, is 
the maximum distance that water travels if the layer 
has a thickness, d. That is, ds=2/. Water saturates 
the soil. Apply load. At the very instant that load 
is applied, as noted before, the water in the soil takes 
all the pressure (see fig. 17-A). That is, at t=0, the 
total pressure, p, is equal to p, at any depth, 2 (see fig. 
17-A). At any other time after loading, p,, is less than 
the total pressure, some of which has been transferred to 
the soil skeleton. The condition, subsequent to apply- 
ing load, is represented by figure 17-B. Atany level, z,a 
part of the load, p,, designated as the “hydrodynamic 
excess,’ is carried by water. The load carried by the 
soil skeleton at depth, 2, is the total pressure minus p,. 

At any time subsequent to f=0, p, is equal to 0 when 
z=0. ‘The water escapes instantaneously at z=0. 

The same condition prevails at the lower boundary, 
that is at2——2/,, 9,=0 1.18 greater than 0; At t=0, 
Pw is a function of z 

For present purposes, assume that the total pressure, 
p, (that on solid and that on water taken together) is 
uniform and constant throughout the depth, 2. 


We now consider equation (41) under the following 
boundary conditions. 


=(Acos Kz+Bsin Kz) 


Pores i Pw=0 when z=0 and ¢>0. 
The sine of 0 is 0 and sin Kz=0. 
Therefore, since p,=0 and B sin Kz=0, it follows 


Ns cK2t | 


(41) 


that A cos Kz must also =0. But cos 0=1 and 
valores = Then A must equal zero. So condition 1 


reduces equation (41) to 


Pw= (Bsin Kz) e-eK* (42) 
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_4p -te 4p tht i St. 
co Pw=—e sing) sobs e sin I, fh 
; (2n-+1)2a%et 
4p cto dita wanes (2n+1)2rz 
Te a sin”—py ~~ ------- (46) 
<a Where n is 0, 1, 2, 3, etc. 
D dz Multiply equation (46) through by dz and integrate 
Ww) between the limits, z=0 and z=2l. We have: 
lo ye al _Onret Poy < 
p 1s) a, {, sin Fdet she J | ain one 
0 bs Tv 0 . 0 Zia 
Guster ae Qn-b1)xe 
dla? cats Ne a 
VI a Gbps lr J, sin 9), Oe 
2lo eee) 2lo 
orf Dp ~dze= — au uw | cosy 
Fal = Yea be 0 0_J0 





-FRACTIONS OF LOAD CARRIED BY SOIL SKELETON 


AND BY WATER. 


FIGURE 18.— 


Condition 2: When z=2/, and p,=0, then sin Ke= 
sin 2Ki,=—0. 


Hence Kz must be an integral multiple of 7. 
z=2I, and Kz=2Kl, for condition Pa NAGh iia 





when Kg=2Klj=a7, 29, 3m 2 in, where n is any 
integer from 1 to reo 
20 on  r 
Since 2Kl,=nz7, K=s>=s7 ar =aT 
4 aT, Sh iD Reeny |: 
— ; on? nxt 3 P 
The expression, e °**, becomes e~ 4h . Substitute 


, : we 2m nT 
for A the successive values, BI,’ BI: Oi,’ etc. in equation 
<b 0 2 lo 


(42), add and obtain: 


mect 4nct 


Ala? “alg? 212 sue, 
Pw=ue sing] tae sin Wl, 
nixrct nre 
ries . [a 
a0 65. el SU ees . Pal4a) 
by 
When t=0, equation (48) becomes: 
2 a 3 
iC, ane : Te 
pay sINgy “1 ae atgp + ds yr +. 
_ N12 
+d, sin 7 (44) 
=" 


This is true for all values of z between 0 and 2/,. To 
3, and integrate all 
terms between the limits, 0 and 2l). Then, as usually 
happens in the case of a Fourier series, all terms on the 
right, except the first, vanish, and since p,—a constant 
=p when t=0, 


2lo 2lo we 
then rf sin; siemens |, sin’57 dz eee (45) 
: Alo 


find a,, multiply through by sin 


std da=0, pee, C,= ap 
Tv 


y= 
3 Nt 


where 7 is odd. 


Substituting the values of a, a2, a3, °° 


*@, IM equa- 
tion (43), there results 














Spiga Smears 

se ll) To? : 

eee) 21, _|, 
Spl, uaa (2Qn-+1) 27] % 
CLs | cos Zion 





n=0o 
_ 16plp , a ee 
1 (2n+1)? 
n=0 


This is the total p, at any time greater than zero. 
2lo : 
The total p, at t=0= i pde=7\2\o=—2kn 
J 
since ~»=p when t=o. 


2lo , 
[ o ,dz=total p, at any subsequent time, f. 
Ai) 








210 e 
| Le pete, the percentage of consolidation at 
Of 

time, tf. 

That is 
16 pl, es Clea Caee oe le 
uponend (e UOT an oe 40 es ) 
i 2 ply 
N e@725N eV 





~9 
or g=1—8/n"(e- +S FS +t s+) 8) 


2 
ee 
where A A 


Figure 18 illustrates diagrammatically the load 
fractions carried by the clay skeleton and the water. 


In figure 19, X = load fraction carried by water, and 
Y = load fraction carried by soil. 


lo 20 

The total | Dy» dz and X+ || p dz=the 
0 0 

total pressure. 


4 Spe Bias aiee: dz Ne 


Sodas yay when ft ree 
“rit _Oate 1 (2n-t)txtet 
n=1—8/r'(e 41? +9e 492 + were Onin 410? ) 
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FIGURE 19.—TRiIANGULAR PRESSURE DISTRIBUTION. THE 
Lower TRIANGLE IS THE MIRROR IMAGE OF THE UPPER ONE. 


DERIVATION OF TIME-CONSOLIDATION RELATIONSHIP FOR TRI- 
ANGULAR PRESSURE DISTRIBUTION, PRESSURE INCREASING 
WITH DEPTH 


For a triangular pressure distribution (see fig. 19) 
the two boundary conditions, 





Condition 1: p,=0 at z=0 


and condition 2: dPw__¢ at 2—(, 


dz 
must satisfy equation (41), 
Dy=e-e (A cos Kz+B sin Ka). 


From condition 1, this equation reduces to equation 
(42) as before and 


pee a asin ele. oe. toy 2 (42) 


For the second condition, at z=Jy, 


Dre __ 9-K2u BK cos Kze=0=—e-«*et BE cos Ky. 


dz 
But 
cos Kl,=0 when Kh="> 
where m is any odd integer. That is, 
T On OT mar 
Kop of he ° Se): hy 


Substitute these values of K successively in (42) and 
add the results together; we thus obtain 


ret 9r2ct 


Sie ee BIC Sa etre 
Po=ae ’sina“+ae 1 sin==+-:- 
21, 21, 
Sater Re ee ee ee Oe (49) 


At t=0, equation (49) becomes 


tee - 3mZ _ mre 
Po Sin ods sin rd -++ +a, sin—-_-. (50) 
21, 2l, 21, 


where m is 1, 3, 5, etc., i. e., any odd, positive integer. 
The coefficient, @», is evaluated in the usual manner. 


At ?=0, Po=t(2)=4= between the limits z=0 and 
0 
z=l,. At t=0, between the limits z= and z=2lh, 


ALR 


Pu=filZ)=271— j 
0 
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Figure 20.—Typrs or Pressure DISTRIBUTION IN A Com- 


PRESSIBLE LAYER. 


Then 


If ("pz .. mre ay Die\« Maz 
an=P|_ Lp sae az+ | (2.7) sin 7? de | 


Where m is any odd, positive integer. 
Integrating, tin =— Pt TOr Hi == loos Gel oa li, Ste. 
T 


and 





On = —2P for M=33 re bt, 15, 19, etc. 


2a? 


By substitution equation (49) becomes 











8p, - TE 8p, an. 8 
Pi 42 . 2 ~Pi 4lp? TZ 
=e 9 ° sin -=-—- 6 Sin eh ws ane ol 
Pu ly Om 2, (51) 
ly m=o a _ mnct lg 
SDpNe Le Aly? . MTZ 

Then dz=—-, » ——>—- e sin —— dz 

Pw 1 m? 215 

0 see 0 

mct Qx*ct (2n-+1)*x°ct 
__16pilo Alo? 16pilo . ~ lo? + 16pilo” at 410? 
ar’ 273° wo (Qn 1)8a* 


Where n=0, 1, 2, 3, etc. 





lo 
At i=, Vi (yee and. | 7 (2) de— Dit 
0 


l 


lo 
Dyd2 
@=1—*n,— 
[i@az 


“= 





2ct 9r2ct 

BD (nr Ge ale 
Sh re 4102 410? : 

oy 7? oF e€ + 


_ 257*ct 
41? 
ym yt... 


— 
LN) 
Or 


Dud 


9 
oe act 
or writing N=—,;>) 
Al, 


oe 1 : 1 : 3 
Go 2(e* 55 on T7955 @25N — ; + Me. ) = (52) 


FORMULAS APPLICABLE TO OTHER TYPES OF PRESSURE 
DISTRIBUTION DEVELOPED 


General relations —Figure 20 illustrates five different 
types of pressure distribution. In this figure the con- 
ditions illustrated are for a compressible stratum over- 
laid with sand and resting on impervious rock. 

The percentages of consolidation for the five types of 
pressure distribution shown in figure 21 are expressed 
as follows: 
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FACTOR,N 
FiGguRE 21.—VALUES OF g AND N For DIFFERENT VALUES OF w. 


Case 1.—Uniform or rectangular pressure distribution. 
=e pS! aN i! —9N 1 —25N 
“a=l1 al e toe Tape “hates 


Case 2.—Triangular pressure distribution, pressure=0 
at upper surface. 


32 1 
n=1- 3 (e"—3 ef oN_, rien .) 


Case 3.—Triangular pressure distribution, pressure =(0) 


attlower surface. 
Ae: 
(4 9 ) 


P 
{eae esate rack 
2 


_ Case 4.—Trapezoidal pressure distribution, pressure 
increasing with depth. 


rates 
Ll oN RS 
pe: 


Case 5.—Trapezoidal pressure distribution, pressure 
decreasing with depth. 


pa—S a (p—pda ot 


ee a 
2 


i= Di 
9 


P p 


If the compressible layer is in contact with permeable 
sand layers above and below, then equation (48) holds 
for all of the five types of pressure distribution shown 
in figure 20 and interpolation formulas are not needed 











in this comparatively simple case for, if g,, the per- 
centage consolidation, is plotted against the time, the 
curves for the different pressure distributions will 
theoretically coincide for the same soil. The ultimate 
amount of consolidation may vary with the type of 
pressure distribution but the percentage of the ultimate 
consolidation attained at the end of a definite period 
of time is theoretically independent of the pressure or 
its manner of distribution if the relatively impervious 
stratum is confined between two relatively permeable 
layers, and the rate at which q increases depends 
solely on the permeability, compressibility, and thick- 
ness of the clay layer. 

The expressions for cases 4 and 5, trapezoidal pres- 
sure distributions, are combinations of equations (48) 
and (52) and are obtained by integration according to 
the procedure followed in deriving equation (52). The 
expressions are, case 4, 


pg +502 
“= Wings (53) 
1 
Dr > 
and case 5, 
pa 
5S ee eee (54) 
a Pa 
eS 


Case 3, triangular pressure distribution, u being 
infinite, is a special case of equation (54). For case 3, 
p=p, and hence from equation (54), 


Is=241— 2 ---------------- (55) 


When the expressions on the right of equations (48) 
and (52) are substituted for g, and q2, equation (55) 
is written 

— ky —N 1 —9N iL —25N 

G—2 “| ¢ +5e Tare alee oges 

peice 32) eee —-9N ny = 25 Nr = ee 4 
Fa bea e 57° TTo5° Pe ea oe) 


Formulas for trapezoidal pressure distributions.—With 
reference to figure 20, both q, and q; approach gq; as a 
limit as p; approaches the limit zero. This is also 
evident from equations (53) and (54). Note also that 


pre is the area of the trapezoid, case 4, when /,=1 and 


p— is the area of the trapezoid, case 5, for =1. 


For the condition, ¢4=q=q:, the numerical value, 
Ng, will be between that for N; and Ng, that is, 


N,=N,+d(N2—N)). 
This is equation (25), Part 1. Then 
(L—q) = (1-@) =(1—q) and log, (1—q) 
=log, (l—q)=log, (l—qu) =log, R. 


Neglecting all terms in e except the first in equations 
(48), (52), and (53), 


log, (1=4,) Sleep ey loo eee (57) 


Tv 
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log, (1—q) STs 5 —M=log. Life (58) 
8 16 
ayy +P 
log. (1—qi) =log, - —iV,=log, RK _~ (59) 
pis 


where log, denotes the natural logarithm. Equations 
(57), (58), and (59) are easily solved for N;, No, and 
N, and the values are substituted in the expression, 


N,—N, 








aN wenn eH ----- (60) 
which is equation (25) solved for J. Thus 
i PACA OS HONIG |b See ace (61) 
4 
log Ee 
Since for case 4, i, err ODE? then by substi- 


tution in equation (61) one obtains 


Putts [Sead 


qo enone = --- (62) 
log = 
which is equation (21), Part 1. 
For the condition, g;,=q;=qs, one may write 
, N;s—Ns3 
fr eee ce: (63) 


Equation (63) is equation (26), Part 1, solved for J’. 
Furthermore, 1—q,=1—q=1—2q+@, from equation 
(55),=1—q;. By dropping all terms in e except the first, 


log. I— eee oe = —N,=log, R as before. __(57) 


and 
log. (1—gs)=log, (l—2q+q2) —Ns3 
16 382 
=log,| 3-2 |-Ni=log. ee eee (64) 
In equation (54), 2¢,—¢q3; may be substituted for qo 
to obtain 
(p— pat as 
ak ee Lee a (65) 
p—5 
Then 
8 te UG: 
o—p)a+0( 5-5) . 
eG 108, ) | IV, (66) 
rs 
ee) 


After substituting the values obtained by solving 
equations (57), (64), and (66) a N;, 3, and Ns; in 


z for p, the expres- 





equation (63) and substituting / 
sion obtained for J’ is 


RUBLE C 








ROADS 19 
ie u(r—2)+2 
Ts SSL Gr—=2) (uF) | Las fee: (67) 
le Lago) 
This is equation (22) of Part 1. 
The expression for J was derived by Gilboy.’ The 
expression for J’ was given by Kimball.'* Neither 


Gilboy nor Kimball indicated the derivations and both 
have given negative signs for equations (62) and (67). 
It would of course be impossible to obtain table 3 of 
Part 1 of this paper with negative values for J and J’. 
Since only the first terms in e were used in their deriva- 
tions, the expressions for J and J’ are of course only 
approximate. 

Figure 21 shows curves obtained by plotting values 
for N on log scale as abscissae against corresponding 
values for g plotted as ordinates. The different curves 
are for 8 particular values assigned to u, these being 
u=0, %, %, 1, 2,3,10,and ©. For intermediate values 
of u, values of g and N may be found by interpolation. 


CHARACTERISTICS OF THEORETICAL TIME-CONSOLIDATION CURVE 
DISCUSSED 


_ Equation (48) is one of a smooth curve having several 
important characteristics, one of which is illustrated in 
figure 22-A and 22-B. F pes sare is the curve 


obtained by taking t= WN, that is, ie 


tion (48). Let q, be the percentage econdeaon at the 
end of a definite time interval, 1, and ?t, the time re- 
quired to produce half this percentage of consolidation, 


qi. i=14 


5 < as unity in equa- 


For q=80 percent, figure 22—-A, 


bole ¢ 


1.6., 


minutes and the time, tf, required to effect 5 Dae 80 per- 


cent consolidation is seen to be 0.3 te from the 
time consolidation curve, figure 22—A. The ratios, 


are plotted against gq, in figure 22—B. Within the 


by 


range g;=0 to q=70 percent, the ratio, - ri ; 1s fairly con- 


stant and is approximately equal to 4. his indicates 
that within this range the time-consolidation curve is 
approximately a parabola. For values of qg greater 


t . 
—, increase. An ordinate 


t, 
value of q,=90 percent on the ye curve, 


than 70 percent, the ratios, 
figure 22—A, corresponds to a value of 7 —==5,27, figure 
22-B. The value, t=2.1 minutes, eens to the 
eer 20 
t ian ad —_ . 


Since the relation between q, and JN, is fixed, the ratio, 


value, 


£5.27 for ae percent is invariant for all values 
1 


assigned to 4 a 5° 

As a rule the time-compression curves for relatively 
impermeable soils have very nearly the same shape as a 
theoretical curve up to a point which, on the theoretical 
curve, has an ordinate of 90 percent. Beyond this 


point the two curves diverge appreciably. Since an 
13 See equation (31) in ‘‘Earths and hela rgs th Proceedings, American Society 
of Civil Engineers, vol. 59, no. 5, May 19 
14 Proceedings, American Society of Civil Engineers, vol. 59, no. 6, August 1933. 
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fo) ordinate of 90 percent on the theoretical curve corre- 
oh sponds to a definite value of pba , lt is easy to com- 
1 
a pute the coefficient of consolidation of a soil by adapting 
# the actual curve to the theoretical one according to the 
ee ay following procedure: 
S 1. Plot the laboratory time-consolidation curve, 
a 60 corresponding to figure 22—A. 
Se ‘ Cea 2. From the time-consolidation curve plot the : ratios 
80 . 4 he 
es corresponding to definite percentages of consolidation. 
100 3. From the f, qm curve find the value of q% corre- 
1 
: t 
TIME: 2 — MINUTES sponding to fae 
A- CONSOLIDATION CURVE 4, From the time-consolidation curve find the value 
of t corresponding to the value of q; for 15,27. Call 
] 
o) this particular value fg. 
5. Substitute the value for ft, thus found in the ex- 
2 
= =~ pression, ne t,=2.1, and solve for c, 
Z 0 
uJ Zeal 0.851)? 
40 eee SS Oe 
= 1ts ty cana Aah) oe (68) 
Ww : : : : 
We sae In testing with the Terzaghi compression device, 
or the sample is confined between two permeable layers 
1 , 
80 and hence b=5 dy. In this case, 
90 
_0.85d,? 
100 CS ie oe (69) 
Daisy 2 This is another method of computing c. [See equa- 
t | tion (8).] 


t 


/ 
B- CURVE SHOWING RATIO = 
l 


FIGURE 22.—RELATIONS BETWEEN PERCENTAGE OF CONSOLI- 
DATION AND TIME, WHEN t=N. 
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